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Abstract. We present the first cosmic shear measurements obtained from the TOOOl release of the Canada-France- 
Hawaii Telescope Legacy Survey. The data set covers three uncorrelated patches (Dl, D3 and D4) of one square 
degree each, observed in u* , g' , r' , i' and z' bands, to a depth of i' — 25.5. The deep, multi-colour observations in 
these fields allow for several data-quality controls. The lensing signal is detected in both r' and i' bands and shows 
similar amplitude and slope in both filters. _B-modes are found to be statistically zero at all scales. Using multi- 
colour information, we derived a photometric redshift for each galaxy and use this to separate the background 
source sample into low-z and high-z subsamples. A stronger shear signal is detected from the high-z subsample 
than from the low-z subsample, as expected from weak lensing tomography. While further work is needed to model 
the effects of errors in the photometric redshifts, this result suggests that it will be possible to obtain constraints 
on the growth of dark matter fluctuations with lensing wide field surveys. The combined Deep and Wide surveys 
give as = 0.89 ±0.06 assuming the Peacock & Dodds non-linear scheme (P&D), and as = 0.86 ±0.05 for the halo 
model and flm = 0.3. We assumed a Cold Dark Matter model with flat geometry and have marginalized over 
the systematics, the Hubble constant and redshift uncertainties. Using data from the Deep survey, the la upper 
bound for wo, the constant equation of state parameter is too < —0.8. 



1. Introduction 

Cosmological weak lensing, also called cosmic shear, can 
be used to probe the dark matter distribution in the 
universe. Weak lensing observations complement other 
probes such as CM B anisotropies (|Spergel et al. 2003|), 
type la supernovae l|Riess et al. 20041 IRiess et al. 19981 
IPerimutter et al. 1999|l . and redshift surveys 
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* Based on observations obtained with 
MegaPrime/Megacam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope 
(CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Science de 
rUnivers of the Centre National de la Recherche Scientifique 
(CNRS) of France, and the University of Hawaii. This work 
is based in part on data products produced at Terapix and 
the Canadian Astronomy Data Centre as part of the Canada- 
France-Hawaii Telescope Legacy Survey, a collaborative 
project of NRC and CNRS. 



l|Lahav fc Suto 2004|l . Weak lensing also has 
the advantage of being free of any assump- 
tion regarding the light versus matter distribu- 
tions HMellier 1999L IBartel mann fc Schneider 200X1 
IVan Waerbeke fc Mellier 20031 [HeFegier 2003| ). 

It has only recently been shown that cos- 
mic shear measurement is technically feasible 
JBacon, Refregier fc Ellis 2000| IKaiser et al. 2000L 

IVan Waerbeke et al. 2000l IWittman et al. 2000|l . 

Unfortunately, the deepest weak lensing survey has 
a sky coverage limited to less than one deg^ and the 
widest to ~10 deg^. Moreover, most surveys were 
performed in one colour only, and even rough redshift 
information was not available. These limitations restricted 
the use of weak lensing as a cosmological probe to a very 
small number of parameters. Early weak lensing surveys 
were primarily focused on the measurements of the nor- 
malization of the dark matter power spectrum, erg, and 
the mass density parameter, flm- The most recent cosmic 
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shear surveys reach a relative accuracy of about 10% 
on a^-^ HMellier 19991 IBartelmann fc Schneider 200ll 
IVan Waerbeke fc Mellier 20031 [Refregier 20"03| ), but the 
uncertainty on other parameters is still fairly large. 

Second generation cosmic shear surveys are now un- 
der way and will provide the community with multi- 
colour data of excellent image quality, over a wide 
field of view. The Canada-France-Hawaii Telescope 
Legacy Survey (CFHTLS)^, using the recently built 
MegaPrime/Megacam wide field camera, belongs to 
this generation. The CFHTLS-Wide survey (the core of 
the CFHTLS cosmic shear survey) will provide a large sky 
coverage of 170 deg.^, and the deep four deg.^ CFHTLS- 
Deep will provide shear information on smaller scales and 
as a function of lookback time, out to higher redshift than 
the CFHTLS-Wide . 

Both surveys will ultimately consist of complete and 
homogeneous panchromatic data in u* ,g' ,r' ,i' ,z' . The 
data were taken between June 1st 2003 and July 22, 
2004 as part of a preliminary survey to provide detailed 
quality assessments and propose technical or operational 
improvements, when necessary. The CFHTLS-Wide has 
the same depth as ViRMOS-Descart, but so far, the 
available data is only in one colour and covers « 20 deg^. 
A cosmic shear analysis with the wide data is performed 
in IHoekstra et al. 20051 The CFHTLS-Deep has been 
observed in all five filters, therefore photometric redshift 
are available, and will be used in this work. Moreover, the 
CFHTLS Dl Deep field is located in the ViRMOS-VLT 
Deep Survey (VVDS) F02 field which has several thou- 
sand of galaxy redshifts (|Le Fevre et al. 2005|l and near 
infrared data (on a tiny area) . A combination of large and 
small scales from the Wide and the Deep data will ulti- 
mately provide an excellent data set to probe the nature 
of dark energy in the universe (|Coora^^^^itere^^99 
IBenabed & Bernardeau 20011 ILinder & Jenkins 2003 



IBenabed fc van Waerbeke 20041 IJarvis et al. 2005|l . In 

this work, we describe the first CFHTLS cosmic shear 
studies based on Deep data, and then combine the Wide 
and Deep data analysis to derive constraints on and 

0-8. 

The organization of the paper is as follows. In section 
2 we introduce the notation and define the statistics we 
use. The data set is described in Section 3. In section 4 
and 5 we present results and residual systematics and we 
discuss them. Conclusions and perspectives are outlined 
in Section 6. 



2. Theoretical background 

The theory of weak lensing has been previously 
been discussed in detail in the literature, in- 
cluding the physical motivations of various ap- 
proximations (e.g. IBartelmann & Schneider 2001). 
Following IHoekstra 2004 Van Waerbeke e t al . 200.2. 



jVan W aerbeke, Mellier fc Hoek stra 2005| we follow the 
notation of Schneider et al. 19981 

We introduce the power spectrum of the convergence 
K as : 
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were Jk (x) is the comoving angular diameter distance out 
to radial distance x{'^)i ^^'^ "-(x) is the redshift distribu- 
tion of the sources. P^d ( fKix) ' 3-dimensional 
mass power spectrum, and k is a 2-dimensional wave vec- 
tor perpendicular to the line-of-sight. 

Cosmic shear can be studied using three different 
2-point statistics, which differ only by their filtering 
schemes. These various statistics have different wavelength 
sensitivities to the power spectrum and therefore the ef- 
fect of systematics on each is different. This enables the 
comparison of multiple cross-checked solutions. Two-point 
statistics are measured as a function of scale 0c, which 
could either be a galaxy pair separation or smoothing win- 
dow radius. The relation between each two-point statistics 
and the power spectrum of the gravitational convergence 
{i.e. the projected dark matter power spectrum) can be 
expressed as follows: 

— Top-hat variance: 
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— Shear correlation function 



(Oe.^^J^ dk kP^{k)Jo{k0,). 
— Aperture mass variance: 

with the aperture mass variance defined as: 



Map(^?c) 



A^e K{e) u{e), 



where U{6) is a compensated filter such as 
9 ^. 0'^\ fl «2 
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Afap can be expressed in terms of the tangential shear 
component inside a circle as follows (fKaiser et al. 19911 
ISchneider 1996|l : 



Afap(ec) 



d^e Q{e), 



(7) 



where the tangential shear component 7t(0) at the posi- 
tion is given by: 
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7t(0) = -Ue (7(0))e-2^^ 



(8) 
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and the function Q{9) is defined as : 

Q{d)^^ J' dO' 9' U{0')~U{e). (9) 

The aperture mass statistic as a tool 
for the cosmic shear analysis has been dis- 
cussed in many papers l|Schneider et al. 19981 

r^en, van Waerbeke fc Mellier 2002} 

Munshi &: doles 2003, Munshi fc Valageas 20051 ). This 
statistic is sensitive to curl-free correlations (E-modes) 
generated by the (scalar) gravitational potential. Curl 
correlations (B-modes) are then easily derived using the 
same statistics, after rotating each galaxy by 45 deg. If the 
only signal present is due to lensing, then the B-modes 
should be zero at all scales. This simple procedure is 
therefore a powerful diagnostic tool to assess systematic 
residuals in cosmic shear signal. 

Unfortunately, the Map statistic is sensitive to the 
smallest accessible angular scales, where cosmic shear sig- 
nal depends on the poorly-known non-linear evolution of 
the dark matter power spectrum. This shortcoming forces 
us to compute E- and B-modes on larger angular scales 
in a different way. For this we use the top-hat shear vari- 
ance and the shear correlation functions. These functions 
are usually derived from the ^_|_ and ^_ shear correlation 
functions: 

C-M = MOhtie + r))-M0M9 + r)), (10) 

where jt and are the tangential and radial projections 
of the shear onto the local frame joining two galaxies sep- 
arated by a distance r. FoUowing lCrittenden et al. 2001al 

we define 

f°° dr' 1"°° dr' 

^'(r)=Ur)+^ ^^_(r')-12r2 / ^^r')- 

(11) 

The E and B shear correlation functions are given by 

kM+m = k(l)_iV). 

A similar relation can be found for the aper- 
ture mass and the top-hat statistics as showed 
in lOittenden et al. 2001 hi lOittenden et al. 2001 hi also 
pointed out and can only be derived up to an inte- 
gration constant which depends on the extrapolated signal 
outside the measurement range. 

Finally, the amplitude of the lensing signal depends 
on the galaxy redshift distribution n{z) (see Eq.QJ). 
As in previous works (see IVan Waerbeke et al. 20021 
|Van Waerbeke, Mellier fc Hoekstra 2005| ), we use the fol- 
lowing redshift distribution: 







) exp 









where a, P and parameters are derived from deep pho- 
tometric redshift catalogues. The lensing signal can be pre- 
dicted for any redshift range using Eq.l^ and Ea. (|13|) . 

3. The Deep CFHTLS TOOOl data set 

The Deep CFHTLS data used in this work consists of 
u* , g', r', i' and z' stacked Megacam images that form 
the first CFHTLS release (hereafter TOOOl). The release 
is composed of stacked images, catalogues and relevant 
meta-data produced from observations in four uncorre- 
lated fields that were carried out at CFHT with the 
Megaprime instrument between June 1st 2003 and July 
22, 2004. Details regarding each field are listed on the 
CFHTLS web pages^. 

Each Megacam image consists of an array of 
9x4 EEV CCD s of 2048 x 4612 pixels each 
HBoulade et al. 2003|l . The pixel scale is 0.186" and the 
camera covers a total field of 1 degree x 1 degree. There 
are two large gaps of 82 arc-second between rows of CCDs. 
In order to produce complete fields, the gaps have been 
filled by organizing observations in a series of exposure se- 
quences with large offsets. This results in an heterogeneous 
pixel illumination at the borders of each CCD. This spa- 
tial fiux variation induces a varying pixel signal-to-noise 
ratio that is taken into account by using pixel weight maps 
together with hand-made masks (see Section 4) to discard 
noisy areas of each field. 

The stacks include only Megacam images with see- 
ing better than 1.0" ^ and airmass below 1.4 have been 
selected. However, because there were fewer u*-band im- 
ages than for the other filters, we relaxed the selection 
criteria for this filter and kept all u* images with seeing 
below 1.4". Only three of the four Deep fields have been 
selected for cosmic shear studies. The D2 Deep field has 
been dropped from our sample because it is significantly 
shallower than the other three fields. 

Data were calibrated and processed at CFHT and the 
Terapix data center. The full TOOOl release is archived 
at CADC^ and available to any CFHTLS registered user. 
A description of the data processing pipeline used to pro- 
duce the deep TOOOl stacks is beyond the scope of the 
paper, but the details can be found on the Terapix web 
pages^. Photometric and astrometric methods and quality 
assessments done on these data are explained in a short ex- 
planatory supplement^. The processing (astrometric and 

^ http : //www . cf ht . hawaii . edu/Sci ence/CFHLS/| 
cfhtlsdeepwidef ields .html 

We use the seeing definition of Terapix as twice the me- 
dian fiux radius of a selection of point sources on each CCD. 
Flux_radius as measured by SExtractor, is the radius of the 
disk that contains 50% of the total flux. For a Gaussian profile 
the SExtractor seeing is almost equal to the Full Width at 
Half Maximum (FWHM). For a typical Megacam PSF, it is 
slightly larger (10%) than the true PSF FWHM. 

http : / / cadcwww.dao .nrc . c a/cfht/cfhtls/J 

^ http : //terapix . iap.fr/ article .php?id_article=382 
® [http : //terapix. iap .fr/article .php?id_article=383 
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photometric calibrations, pixel re sampling, image warp- 
ing and stacking, catalogue production) uses the current 
first generation Terapix software tools and closely follows 
the one used for the ViRMOS-Descart survey that is de- 
scribed in lMcCracken et al. 20081 We refer to this paper, 
and to the Terapix and CFHT'^ web pages for further 
details. 

The accuracy of the photometric calibrations can be 
estimated from the stellar colour-colour plots and the 
galaxy counts in all bands given on the Terapix TOOOl 
pages and is also discussed in the more detailed stellar 
analysis done by Schulteiss et al. (in preparation). In all 
bands, the cumulative internal and systematic photomet- 
ric errors are 0.05 mag up to AB=22.5, and never larger 
than 0.1 to the 80% completeness hmit (« AB=25.5). This 
uncertainty is sufficient for the cosmic shear studies on this 
paper. 

Table H summarizes the TOOOl stacks used in this 
work. The completeness limits have been computed by 
adding randomly simulated stars (Moffat profiles) inside 
a 2000x2000 area of each Deep field and by running the 
detection and photometry again, using the MAG_AUTO mag- 
nitude of SExtractor* software (Bertin & Arnouts 1996(1 . 
The completeness was also checked using galaxy counts ^ . 

The D4-z' data have not been released by Terapix 
because the astrometric internal accuracy was below the 
scientific requirements. The large rms error found in 
the HA-z' data is not been fully understood, but it re- 
sults in a large number of galaxy mismatches during the 
u*, g', r', i', z' catalogue cross-identification. This hampers 
reliable panchromatic studies for many galaxies detected 



are then used to derive orientations and raw ellipticities 
of galaxies from the weighted second moments of the 
galaxy light distribution. In order to minimize the noise 
contribution each moment is filtered using a Gaussian fil- 
ter W{9) of size r/i! 



d?ew{e) 0,9,(9) f{9) 



where f{9) is the surface brightness. 
The raw ellipticity is given by: 



(14) 
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where Tr{I) represents the trace of the matrix /. We use 
the KSB method to get an unbiased estimator of the shear 
7. This method has been tested by several teams and it 
has been demonstrated that it provides robust and re- 
liable shear measurements from ground based data (see 
the comprehensive critical investigation of KSB and other 
techniques by [Heymans et al. 2005b[ and also references 
therein). 

4.1.1. PSF correction: the principle 

Let us assume the shear-free intrinsic ellipticity of a galaxy 
is e''. On the detector, its shape is eventually modi- 
fied by the distortions produced by gravitational lensing 
effects and systematics that increase the smearing and 
the anisotropic component of the PSF (atmosphere, op- 
tical aberrations). Assuming these distortions are small, 



in this field. A further investigation reveals that the D4-i' iKaiser et al. 1995, demonstrated the observed ellipticity, 



astrometric solution is also slightly off, while data quality 
in other filters are excellent. Although it has no impact 
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photometric studies, a quick weak 
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lensing analysis of the D4-i' field shows it has more sys- 
tematic residuals than Dl-j' and D3-i'. In contrast, the 
three Deep r' band data have similar quality and do not 
show systematics residual differences. We therefore used 
the r' band as the reference data set for all comparison 
between the fields, and only use the deep D1/D3 i' band 
data for colour comparisons, when needed. 



where q is the anisotropic component of the PSF and 7 
is the gravitational shear. P'^^ and P'*™ are called the 
shear and the smear polarisability. Their values depend 
on the galaxy surface brightness and on the filter proper- 
ties W{9). q can be derived directly from the data, by 
measuring the ellipticity of stars in each field, e^,, such as: 



4. Detection of the shear signal 

4.1. Galaxy shape parameters 

Catalogues and shape measurements of galaxies are pro- 
duced using the IMCAT software UKaiser et al. 19951 here- 
after KSB). For each object the centroid position and 
the half-light radius r/j are measured. These parameters 



http : //www . cf ht . hawaii . edu/Science/CFHTLS-DATA/ 1 
dataprocessing . html 

* http : //terapix . iap . f r/rubrique . php?i d_rubrique=91| 

® see http ://clix. iap. fr/TOOOl/Plots/j 
CFHTLS_D_i_galcouiit_T0001 .png 
and http : //clix . iap . ir/T0001 /Plots/| 
CFHTLS_D_r_galcount_T0001 .png 
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The shear polarisabiHty is however ahered by the 
isotropic smearing component of the PSF. It results in 
a modification of the shear polarisabihty 



JDsh 

p-y psh -A- psm 

r>sm 



(18) 



where is called pre-seeing shear polarisability and 
psh/sm j.gfgj.g stars ( |Luppino fc Kaiser 1997| ) . Provided 
the assumption (eo) = is valid, an unbiased estimator of 
the shear 7 is given by: 



7=( 



P""q)). 



(19) 
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Table 1. Summary table of TOOOl Dl, D3 and D4 deep stacks used in this work. Magnitudes are instrumen- 
tal AB. Details on magnitude, aperture, seeing and completeness definitions are given in the explanatory page 
jhttp : //terapix . iap . f r/article . php?id^rticle=383| 





Dl 


D3 


D4 


RA (J2000) 


02:25:59 


14:19:27 


22:15:31 


DEC (J2000) 


-04:29:40 


+52:40:56 


-17:43:56 


Effective FOV (deg^) 


0.80 


0.77 


0.77 


Exp. time u* (s) 


10560 


4620 


16680 


Median seeing u* (arc-sec.) 


1.15 


0.88 


1.05 


Completeness u* 50% (mag.) 


26.4 


26.0 


26.2 


Exp. time g' 


7515 


8010 


11250 


Median seeing g' (arc-sec.) 


0.98 


0.95 


0.99 


Completeness g' 50% (mag.) 


26.4 


26.5 


26.2 


Exp. time r' 


17280 


20820 


26400 


Median seeing r' (arc-sec.) 


0.87 


0.93 


0.85 


Completeness r' 50% (mag.) 


26.1 


26.4 


25.9 


Exp. time i' 


52000 


59640 


58800 


Median seeing i' (arc-sec.) 


0.88 


0.92 


0.88 


Completeness i' 50% (mag.) 


26.1 


26.2 


25.8 


Exp. time z' 


12240 


15120 




Median seeing z' (arc-sec.) 


0.86 


0.85 




Completeness z' 50% (mag.) 


24.5 


24.6 





4.1.2. Object selection 

Prior to cosmic shear analysis, all CFHTLS images are 
checked by eye and masks are drawn by hand. These 
masks are designed to avoid elongated defects, like sat- 
urated stars, as well as large foreground galaxies with ex- 
tended bright halo that may contaminate the shape of un- 
derlying faint galaxies (see IVan Waerbeke et al. 200 1" for 
details). We should emphasize that masks are only drawn 
using criteria (2; = galaxies, bright stars, CCD defects) 
that are not correlated with the lensing signal. In addition, 
we used the weight map images produced by Terapix for 
each stack to reject all pixels with a relative weight ampli- 
tude less than of 80%. This rejection step reduces signifi- 
cant spatial variation of the detection threshold and keeps 
the averaged redshift distribution of lensed galaxies sta- 
ble over the field. The rejection scheme removes the CCD 
boundaries from all of the fields, and is essentially equiva- 
lent to singling out each CCD region, as was done earlier in 
ViRMOS-Descart survey. The gain in homogeneity is how- 
ever preserved at the expense of the sky coverage. About 
30% of the initial area is lost after the masking process. 

Stars needed for the PSF correction are selected 
along the stellar locus of the magnitude/size diagram 
IjFahlman et al. 1994|, from the region where stars are 
about one magnitude fainter than the saturation level 
and where they cannot be confused with faint galaxies. 
The Pf™ and P^'' values are derived at all Megacam im- 
age positions from a PSF mapping that samples the PSF 
smearing and PSF anisotropy at the position of each star, 
and by interpolating their values between the stars. This 
operation is done on each CCD separately, as suggested 
bv lHoekstra 20041 The PSF is mapped using a composite 



model of a second order polynomial and a rational func- 
tion, p°'{x,y): 

P"{x.y)^ (20) 
ao + aix + a2y + a-^x^ + a^xy + a^y^ + c{x, y) 

where c(x, y) is the rational function chosen as: 

c{x,y)^ (21) 
bp -I- hix + h2y + b^x"^ + bjxy + b^y"^ + bgy^ + 67^/"^ 
l + bsx + bgy 

The second order polynomial terms models the smooth 
low frequency PSF component, while the rational func- 
tion provides a model for the high frequency PSF terms 
(Hockstr a 2004|l . 

The correction is made in two steps. First, the 
coefficients of the rational function are determined. Since 
the CFHTLS Deep fields are much deeper than the 
RCS IjHoekstra et al. 2002jl and the ViRMOS-Descart 
(|Va,n Waerbeke et aJ. 20001 IVa,Ti Wa,erheke et al. 2001 |l 
surveys, the density of selected stars is higher and we do 
not need to map the PSF using external stellar fields. 
Each field has about 100 stars per CCD, so the high 
frequency PSF terms can be reasonably well sampled 
down to 0.5 arc-minute, and all coefficients of the rational 
function can be constrained with sufficient accuracy. In a 
second step, the polynomial terms are determined. 

We also compared the rational function solution 
against the second order polynomial interpolation. We 
found the results are not very different from our com- 
posite model, although the rational function improves the 
quality and stability of the PSF mapping. 
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Magnitude 

Fig. 1. Magnitude distribution of galaxies 

in the three fields in the r band. More de- 
tailed galaxy count plots, for each filter and for 
each Deep field, are available on the web at 
|http : // terapi x . iap . f r/article ■php?id^rticle=382| 



Once ellipticity is corrected we keep in the sam- 
ple all objects with angular size larger than the see- 
ing disk and smaller than two arc-seconds. Following 
IVan Waerbeke et al. 2000l we automatically reject one 
galaxy in every close pairs with angular separation less 
than 12 arcsec in order to avoid contamination of elliptic- 
ity measurements by overlapping isophotes of neighboring 
galaxies. 

The magnitude distribution of the final object cata- 
logue is shown in Fig^ The limiting magnitude corre- 
sponding to a 80% completeness hmit is r^^ = 25.5. 
Bright objects with magnitude smaller than 21.5 and faint 
objects with magnitude larger than 25.5 are also removed 
from the galaxy sample. The final galaxy number density 
of the cosmic shear catalogue is about 20 / arcmiri^ . 

As proposed bv lErben et al. 20011 we assign an ellip- 
ticity dispersion ag to each object corresponding to the 
ellipticity dispersion in a box containing its 20 nearest 
neighbors in the {magnitude, size) space. Weighted 2- 
point statistics are computed assigning to each galaxy a 
weight given by -I- af) where CTe is the ellipticity 

dispersion of the unlensed galaxies. A different noise esti- 
mation HHoekstra et al. 2000|l gives similar results. 



5. Residual systematics 

5.1. Quality of the PSF correction 

A visual inspection of the MegaPrime PSF (FigEJ shows 
that the PSF anisotropy has significant variation over the 



field and may also be very large at the boundaries^". The 
PSF correction is therefore a critical step and its reliability 
demands careful verifications. In addition to the usual B- 
mode analysis shown in the next section, in this section 
we carry out several analyses of the systematics. 

The quality of the PSF correction and its homogene- 
ity over the Megacam field camera can be assessed by 
comparing the mean star ellipticity before and after PSF 
correction (FigOl). The average stellar ellipticities (et) and 
(e,-) are plotted as a function of the radial distance from 
the center of the field, r. It is interesting to note that 
the radial ellipticity component degrades much more and 
much faster than the tangential ellipticity. However, the 
PSF correction done by the PSF mapping is very good, 
for each Deep field. After correction, the dispersion of star 
ellipticities is about 2 x 10""^ at any point of the camera. 
There is no significant change in the residual error as func- 
tion of position. The small increase in the fluctuation of 
star shapes at very small distances is due to higher Poisson 
noise: each radial bin has the same width, so the inner- 
most circle encompasses the smallest area and contains 
fewer stars than the others. 

Finally, we checked the residual amplitude of the shape 
correlation function between corrected stars. We found is 
to be two order of magnitude smaller than the expected 
lensing signal at all scales probed by this work (Fig.^. 

The tests discussed above only guarantee that the PSF 
correction is excellent in the neighborhood of selected stars 
or on angular scales larger than, or close to, the mean an- 
gular distance between stars. In regions where no stars 
were selected or on small scales, the local PSF correc- 
tion residuals may be larger than the average. A use- 
ful test of systematic residuals on small scales has been 
proposed bv'Baco n et al. 2003l and'Hey mans et al. 2005b| 
Assuming the PSF model derived from stars and applied 
to galaxies is unable to remove all systematic contribu- 
tions, the star-galaxy cross correlation will be non-zero 
and may vary as function of angular scale. If the residual 
is small. iBacon et al. 2003l showed the systematic residual 
can be expressed as follows: 

(e. e.) ' ^ ' 

where egai is the corrected galaxy ellipticity and e* is the 
uncorrected star ellipticity. We use the ^sys to compute 
the contribution of systematics for both top-hat and com- 
pensated filter. Fig. |31 shows they are consistent with zero 
at all scales between 0.5 arc-minute to 30 arc-minutes. 
This confirms that residual systematics are negligible in 
the Megacam Deep fields. 

5.2. Independent analysis of r' and i' data 

The robustness of cosmic shear signal can also be assessed 
by comparing results obtained using different filters for 

This strong PSF anisotropy has been considerably reduced 
by the CFHT staff, after the TOOOl release. It should no longer 
be a critical issue for next releases. 
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Fig. 2. The mean ellipticity of the stars shows the behavior of the PSF anisotropy for the three fields Dl, D3 and D4 
. For all of the three fields, < e > is few percent in the central part and it becomes about 10% in the corners (see also 
fig©. 
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Fig. 3. Bottom panels show the mean tangential (left 
panel) and radial (right panel) uncorrected stellar elliptic- 
ity as a function of the distance to the center of the camera 
for DlrTOOl (filled triangles ), DSrTOOl (filled squares) 
and D4rT001 (filled circles) fields. Top panels show the 
same quantities after PSF correction. 



the same galaxy sample. Because gravitational lensing is 
achromatic, we expect the shape and amplitude of cos- 
mic shear to be identical for data taken in different filters. 
Any significant difference between two bands provides a 
diagnostic of the PSF corrections. A first attempt at com- 
paring shear measurements in different filters was made 
bv lKaiser et al. 20001 using the CFHT12K camera. The / 
and V bands showed significantly different signals that 
were inconsistent with the change in redshift distribution 
between the two filters. 
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Fig. 4. E-modes (red filled circles ) and B-modes (black 
open circles) top-hat two point statistics of corrected stars 
show the smallness of residual PSF systematics. 



The Deep photometry provides a sample of the same 
galaxies detected in different filters, so the signal is ex- 
pected to be the same. However, these filters have differ- 
ent depths, and the shallowest colours do not have enough 
galaxies to allow a comparison of the signal between all 
colours using the same galaxies. This limitation affects 
mainly the u* and g' bands. Furthermore, these bands are 
more sensitive to atmospheric dispersion than other filters. 
We expect their PSF anisotropy to be larger than for r', i' 
and z' bands and its correction may also depend more on 
the relative differences between the averaged spectral en- 
ergy distributions of stars used for the PSF calibration and 
of galaxies. Hence, u* and g' are not well suited for weak 
lensing analysis. The comparison between the r' and the 
z' bands doesn't give many informations because of the 
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The r' and i' bands results are remarkably similar, both 
in shape and amplitude, they agree to within la at all 
scales. 
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Fig. 6. Top-hat variance of the E-modes in the i'-band 
(red filled circles) and in the r'-band (red filled triangles) 
for the same data set. Top-hat variance B-modes in the 
i'-band (black open circles) and in r'-band (black open 
triangles). 



Fig. 5. Two-point signal statistics(red filled circles) :top- 
hat (bottom) and Map (middle) correlation function (top) , 
compared with residual systematics (black open circles). 
Signal error bars are statistical ones. 



limited size of the matched sample. We therefore decided 
only to focus on the comparison between r' and i' bands. 

We computed the two-point statistics using the same 
objects in i'-band and r'-band in the Dl and D3 fields only. 
As reported before, D4 was discarded from this study be- 
cause it shows higher systematic residuals in i' bands than 
the two other fields. It is worth noting that both r' and i' 
band images have been processed (fiat fielding, astromet- 
ric and photometric calibrations, image selection, image 
stacking) in a totally independent way. The only correla- 
tions between the two samples are the software tools and 
the pipeline scheme used at Terapix. 

The r' and i' ellipticity catalogues have been computed 
and PSF-corrected independently, starting from the r' and 
i' TOOOl stacked images. The galaxy cross-identification 
is done at the very end of the processing to compare the 
results. Fig. IHl shows the comparison of the E- and B- 
modes of the top-hat shear variance for both the i' and r' 
data sets. The error bars are estimated as the quadrature 
sum of the statistical and the systematic error S^sys defined 
by Ea. f^ . The amplitude of the latter is bigger in the 
r' band as shown by the residual B-modes in this filter. 



6. Characterization of the shear signal 

6.1. Two point statistics 

The ellipticity correlation functions £,+ {r) and £,-{r) are 
measured from the weighted mean of all pairs with an- 
gular separation r. The correlation function is computed 
using equations and H12|) . The Map and the top- hat 
statistics are also computed as a function of the correlation 
functions ^+(r) and ^_(r) following ISchneider et al. 20021 
and ICrittenden et al. 2001 bl FigEI shows the two-point 
statistics for the three deep fields Dl, D3 and D4. Error 
bars including statistical noise and cosmic variance are 
computed from the S,+ {r) and S.-(r) as described in 
(Schneid er et al. 20021 

The cosmic variance contribution is computed using 
the CFHTLS TOOOl Deep survey properties: an effective 
density (after masking) of 20 gal / arcmiri^ , an effective 
area of 2.1 deg^, and an ellipticity dispersion per elliptic- 
ity component of 0.3 (the latter was measured from the 
corrected ellipticity). However, the error calculation de- 
scribed in ISchneider et al. 20021 is only valid for a single 
connected field with a number density of n equally-sized 
galaxies. We therefore replace the statistical error compo- 
nent by the Poisson noise measured from the data, using 
the weights (computed as described above) and positions 
of each galaxy. For the top-hat variance and the corre- 
lation function, the free integration constant is chosen so 
that the B-modes on scales between 15 and 25 arc-minutes 
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vanish. Although its amphtude is meaningless, one can see 
that the B-mode is flat and stable over that range of an- 
gular scales. 

In contrast, the Map statistic does not have an undeter- 
mined integration constant (as explained previously), so 
the B-mode amplitude is a physical property. Fig|3 shows 
the presence of B-modes. Note that the Map filter for a 
given size is mostly sensitive to scales around ~ 6*75. 
This explains why the other two-point statistics do not 
show B-modes at the same scales. The B-mode at such 
small scales may result from intrinsic alignment of galax- 
ies dKing fc Schneider 2002| [Heymans fc Heavens 20031 ) 
or from the correlation between intrinsic ellipticity and 
shear ( Hirata fc Seljak 2004| ). If these systematics are real, 
we expect to correct them in future work by using the pho- 
tometric redshifts. A further investigation confirms that 
the B-modes come from weak objects (i.e. 25.0 < < 
25.5), and that a magnitude cut that rejects objects with 
magnitude fainter than 25.0 gives zero B-modes at all the 
scales, even for the Map statistic. However, we keep these 
objects in our catalogues because a deep sample will be 
necessary to study the evolution of signal with redshift. 
In addition, the presence of B-modes at small scales will 
taken into account when we estimate cosmological param- 
eters. 



6.2. Evolution of signal with redshift 

The cosmological nature of the two-point statis- 
tical signal can be established by comparing its 
amplitude as function of source redshifts with the- 
oretical expectations of the gravitational instabil- 
ity paradigm and the gravitational lensing the- 
ory IjBernardeau et al. 19971 |Jain fc Seljak 1997| ). 
To first order, the signal should increase as zj^^ 
HBernardeau et al. 19971 |Jain fc Seljak 1997| so even a 
rough separation of galaxies into low- and high-redshift 
populations should split the cosmological lensing signal 
accordingly. 

The CFHTLS TOOOl data sets are well suited for this 
analysis. The observations can be used to sample the high 
redshift universe up to z ~ 1. There are enough of galax- 
ies to divide into two subsets based on their estimated 
photometric redshifts. 

Photometric redshifts were measured 
using the hyper-z public software^^ 
HBolzonella, Miralles fc Pello 2000| ). hyper-z uses the 
multi-band photometric data of a galaxy to derive its 
most likely redshift and spectral energy distribution 
(SED) based on the Bruzual & Chariot evolution models 
IjBruzual fc Chariot 1993}) . 

We used the Dl and D3 u*, 5', r', i' and z' images 
and the D4 u*, 5', r', i' images (the D4-z' stacked im- 
age is missing in TOOOl). Photometric catalogues were 
produced by the SExtractor software. All galaxies were 
first detected in the r' band reference image. Magnitude 
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Fig. 7. Two-point statistics for all the three fields com- 
bined. Red filled circles show E-modes, black open circles 
show B-modes. E-mode error bars include the statistical 
error and the cosmic variance contribution, while B-modes 
are affected only by statistical error. 



and colours of galaxies are then computed using the 
r'— center positions and inside an aperture scaled accord- 
ing to the size of each galaxy in r'-band. The minimiza- 
tion was performed assuming magnitude errors derived 
from SExtractor, which range between Amag — 0.03 and 
Amag = 0.1 in all bands. 

Fig. IHl shows the photometric redshift distribution of 
the galaxies in Dl field down to i' — 24.0. This subsam- 
ple can be compared with the VVDS spectroscopic red- 
shift distribution obtained from 11000 spectra in the same 
region ULe Fevre et al. 2005). There are no apparent dis- 
crepancies that would make the separation into photo- 
metric low- and high-redshift galaxies unreliable. Beyond 
i' = 24.0, large spectroscopic redshift samples are not yet 
available, but we don't have any reason to believe that our 
photometric redshift accuracies will degrade significantly 
for the i' < 24.0 sample. 

The cosmic shear catalogue can therefore be split 
into two samples with equal numbers of galaxies at high 
and low redshifts with reasonable confidence and can be 
compared with cosmological predictions. Poisson noise is 
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Fig. 8. Density of galaxies of Dl field as a function of 
photometric redshift. 



Table 2. Mean photometric redshift in magnitude bins. 



magnitude bin 



18.5 < iAB 
19.0 < lAB 
19.5 < lAB 

20.0 < iAB 
20.5 < iAB 

21.0 < iAB 
21.5 < IAB 

22.0 < IAB 
22.5 < iAB 
23.0 < iAB 
23.5 < iAB 



< 24 

< 24 

< 24 

< 24 

< 24 

< 24 

< 24 

< 24 

< 24 

< 24 

< 24 



mean redshift 



< z >- 

< z>-- 

< z>-- 

< z>-- 

< z >- 

< z >- 

< z >- 

< z>-- 

< z>- 

< z >- 

< z >- 



0.850 
0.853 
0.858 
0.865 
0.876 
0.892 
0.913 
0.942 
0.981 
1.035 
1.100 



Fig. 9. Top-hat variance for "high-z" subsample (red 
filled triangles) and "low-z" subsample" (red filled circles). 
B-modes for the two subsamples are also shown. Error bars 
include statistical noise and cosmic variance. The data 
are compared with theoretical fiducial model {Qm — 0.3, 
Ha = 0.7, as = 0.88 and h = 0.7 ) and source distri- 
bution n{z) modeled by Eq^J with a — 1.98, f3 = 0.66 
and Zs = 0.0981. The low-z source selection is simulated 
using n{z) between 0.3 < z < 1.2 and zero otherwise (bot- 
tom line). Likewise, the red filled triangles and the black 
open triangles are the E- and B-modes of the "high-z" 
sample. The data are compared with the same theoretical 
model with a high-z source selection simulated using n{z) 
between z > 0.8 and zero otherwise. Shaded areas show 
models within Zg — 0.0981 j^Q gn that represent the la er- 
ror region on Zg as derived from the likelihood parameter 
estimation. 



therefore similar in the two subsamples, but photomet- 
ric redshift errors are expected to be larger in the high- 
redshift tail. 

Fig. El shows the top-hat shear variance measured 
for the two populations. The low-z sample ranges in 
0.3 ^ z ^ I., while the high-z galaxies have z ^ 1.0. Error 
bars include Poisson noise and cosmic variance (see Sect. 
4.2). The difference between the two samples demonstrates 
the cosmological nature of the signal. An indicative com- 
parison of signals with theoretical predictions is also plot- 
ted. 

The relative lensing amplitude for the two source 
galaxy populations is less sensitive to cosmic variance fluc- 
tuations, and agrees with the predictions. 

Contamination by galaxies with incorrect photometric 
redshifts is likely important, in particular for the faintest 
galaxies and the high-z tail ( further informations about 
degeneracy of photometric redshifts in the case of missing 
infrared bands can be found on hyper-z user's guide). 



In spite of potential contamination by incorrect pho- 
tometric redshifts, the cosmological imprint of large-scale 
structure detected in the Deep CFHTLS data shows that 
MegaPrime is suitable for cosmic shear studies. Fig. 
El also demonstrates that the CFHTLS Deep survey has 
promising potential for tomographic studies that explore 
the evolution of the dark matter power spectrum with 
look-back time. 

7. Parameter estimation 

7.1. Derivation of the likelihood function 

In this Section, we describe the estimation of cosmological 
parameters. In a subsequent paper, we will perform a com- 
plete parameter estimate, combining weak lensing with 
other cosmological probes. Therefore, here we limit the 
analysis to the normalization of the mass power spectrum 
((Ts) and matter density (flm) measurements. The shape 
parameter P is given by the Cold Dark Matter paradigm 
P = f7„i h, where h is the reduced Hubble constant. We 
allow the characteristic redshift of the source distribution 
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to vary around the best fit that will be described in the 
next sub-section. 

To measure cosmological parameters, we adopt a 
maximum-likelihood method. Let di be the input data vec- 
tor (i.e. the top-hat shear variance as a function of scale 
9i), and mi(r2,„, erg, n(z)) the prediction, fmiction of the 
parameters to be estimated. The likelihood function of the 
data is then: 

^ " (27r)"'|C|V2 ['^'^^ ~ - ^^f] ' (23) 

where n = 16 is the number of angular scale bins and C 
is the 16 X 16 covariance matrix of the top-hat shear, 

Cy = ((d, -m,f K -m,)), (24) 

and C can be decomposed as C — C„ + C^, where C„ is 
the statistical noise and C<; the cosmic variance covariance 
matrix. 

As discussed above, the matrix is computed accord- 
ing to lSchneider et al. 20021 assuming an effective survey 
area of the CFHTLS Deep fields: 2.1 square degrees, a 
number density of galaxies Ugai = 20 /arcmin^, and an 
intrinsic ellipticity dispersion of — 0.3 per component. 

The cosmic variance is computed assuming Gaussian 
statistics. While this assumption becomes inappropri- 
ate on small angular scales, errors on such scales 
are dominated by the statistical noise contribution, so 
the Gaussian approximation remains an excellent one 
IjVan Waerbeke et al. 2002)l . The covariance matrix com- 
ponents are derived for a fiducial cosmological model cor- 
responding to the best fit of WMAP data proposed by 
ISpergel et al. 2003| VIm = 0.3, 17a = 0.7, erg = 0.88, 
r = 0.21 (the reduced Hubble constant is h^O.7). The B- 
mode is calibrated by marginalizing around B — within 
the 1(7 interval. 

7.2. Parameter estimation 

The source redshift distribution is calibrated us- 
ing the Hubble Deep Field (HDF) catalogues 
IjFernandez-Soto et al. 1999jl . which provide a more 
accurate estimate of redshift in absence of infrared data 
in CFHTLS fields. It turns out that the F606 filter of 
WFPC2 is a good match to the Megacam r' filter 
within our 1 a magnitude error. We select all galaxies 
with 21.5 < r' < 25.5. The Hubble Deep Fields provide a 
sample at high redshifts that overlaps with the redshift 
range expected for the CFHTLS Deep fields. 

We use the source redshift distribution model of 
Eq.l2 and perform a fit, allowing the parameter Zg 
to vary. We then identify the ilcr and ±2cr uncertain- 
ties, which we marginalized over in the cosmological pa- 
rameter estimation. We find a = 1.9833, f3 — 0.6651, 
z, = 0-098ltomltom6- Figure CUl shows the unnormal- 
ized weight in magnitude slices in the Deep catalogues. 
The effect of down-weighting faint galaxies is taken into 
account in the source redshift estimation. Fieiire lTTl shows 



the best fit model and the underlying photometric red- 
shifts from the Hubble Deep Fields (solid line) . Error bars 
are Poisson errors. The dashed-dotted line on Figure ITTI 
shows the redshift distribution one would have if we ignore 
the weighting. The best fit redshift distribution model has 
a mean source redshift of « 1.01, nearly 0.2 higher in z 
than the Wide survey (Hoekstra et al. 2005). 

The constraints on flm and erg are obtained after 
marginalization of the reduced Hubble constant h £ 
[0.6, 0.8] and over the ± 2(t limits of the source red- 
shift parameter Zg . The resulting constraints in the flm-o's 
plane are given in Figure ^1 This figure shows that the 
CFHTLS Deep field gives constraints as good as previ- 
ous Icnsing measurements, despite its small field of view. 
This is the consequence of the larger fraction of high red- 
shift galaxies, which are more strongly lensed. Using the 
Peacock & Dodds (1996) non-linear scheme, we obtain 
erg = 0.94 ± 0.15 ± 0.20 (± Icr ± 2a) for f7„ = 0.3. Error 
bars are the one and two a errors respectively. The Smith 
et al. (2003) halo model gives dg = 0.90±0.14±0.20, which 
agrees with previous normalization measurements. The 
similarity between the result obtained using the Peacock 
& Dodds and that using the halo fitted model is not sur- 
prising. Indeed, on scales 1', which dominate our signal, 
the difference between the two models of power spectrum 
is ^ 5%. On smaller scales, we would expect an increasing 
discrepancy between these different ways to estimate erg. 

We then measure erg by combining these constraints 
with those obtained on the CFHTLS Wide survey (see 
Hoekstra et al. 2005 for the details) . The result of this joint 
analysis is shown in Fig ure 1131 and remarkably, the i^m- 
(Tg degeneracy is partially broken. This is the consequence 
of measuring the large and small scales simultaneously, 
as shown in [Jain fc Seljak 19'97| For = 0.3, we get 
erg = 0.89 ±0.06 ±0.12 using Peacock & Dodds (1996) for 
the non- linear scheme and cg — 0.86 ± 0.05 ±0.11 for the 
halo model (Smith et al. 2003). 

The power spectrum normalization is in very 
good agreement with results from medium-redshift 
and low-source-redshift weak lensing surveys 
(Hoekstra et al. 20021 IVan WaerbekT fc Mellier 20^. 
It is remarkable that the parameters of the redshift 
distribution, which have been estimated from a different 
survey, are such that the normalization erg lies within the 
errors of previous measurements. This is strong evidence 
that deep, medium and shallow lensing surveys are in 
cosmological agreement, hence reinforcing the ability of 
cosmic shear to probe the mass distribution at different 
redshifts and different scales. 

Weak lensing can also be used to constrain dark en- 
ergy. Figure [T^ show the upper limit on wq, the constant 
equation of state parameter derived from the Deep data 
only. Here we used only the Peacock & Dodds non-linear 
prescription (a detailed discussion on non-linear power 
spectrum correction in the context of Dark Energy can 
be found in Jlockstra et al. 2005 , which also includes a 
joint analysis of the Wide and Deep data). We obtain 
Wq < —0.8 at ler, and the contours show that this re- 
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suit is independent of . This is particularly interesting 
because lensing combined with either cosmic microwave 
background (Jarvis et al. 2005) or supernovae will provide 
a strong constraint on the dark energy equation of state. 



23 24 
mag r 



Fig. 10. Plot of mean weight per galaxy as function of 



magnitude 21.5 < r^^ < 25.5. 
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Fig. 12. fl„i and as constraints with the Deep data only. 
The contours show 0.68, 0.95 and 0.999 confidence regions. 
Errors include statistical, covariance and residual system- 
atic contributions. The models are pure Cold Dark Matter 
fit to the data, marginalized over the redshift distribution 
(see Section 7.2 for the details). 
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Fig. 13. Same as Figure [HI combined with the CFHTLS 
Wide data (Hoekstra et al. 2005). For ^Im = 0.3 we have 
(Ts — 0.86 ± 0.05 at la (see Section 7.2 for details of the 
error calculation). 



Fig. 11. The histogram shows the photometric redshift 

distribution of 21.5 < r'^g < 25.5 galaxies of Hubble Deep g_ Summary and conclusion 



Field North and South used in this work. The central solid 
line is the best fit model. The solid line histogram is that 
magnitude weighted redshift distribution. The dashed-dot 
histogram shows the redshift distribution if the galaxies 
were not magnitude weighted. 



This paper describes the first cosmic shear studies of 
CFHTLS Deep data using the TOOOl CFHTLS release. 
It uses data collected in u* , g' , r' , i' and z' with 
Megaprime/Megacam over the first year of the survey. 
Only between 1% and 15% of the Deep data are therefore 
in hand depending on the field and on the filter, and so 
the survey is still 2-3 magnitude below the final goal. 
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Fig. 14. Dark energy constraints from the Deep data only. 
Hidden parameters are marginalized using a flat prior over 
as e [0.7,1.0], h G [0.6,0.8] and within the ±2cr bound- 
aries of the redshift parameter Zs (see Section 7.2). 



The TOOOl data have been used to assess the capabil- 
ities of Megaprime/Megacam and to clarify the poten- 
tial and the science drivers of the CFHTLS Deep survey 
for weak Icnsing studies. 

The correction for PSF anisotropy works very well, 
showing that residual systematics are almost zero at 
all scales probed by a Megacam field. This is con- 
firmed by the star-galaxy cross-correlation analysis. This 
also demonstrates that the CFHT-Elixir-Terapix cal- 
ibration/reduction pipelines can deliver co- added images 
which have the required lensing quality. However, the pres- 
ence of B-modes by weak objects at small scales should be 
further investigated. 

The cosmic shear signal has been detected in the r'- 
band. Its consistency and achromaticity has been checked 
by independent r'- and i'- analysis of the same data sets. 
We have presented results for three standard two-point 
shear statistics. 

Thanks to the depth of the CFHTLS Deep sample, 
and using the photometric redshifts derived from the u*, 
g' , r' , i' and z' images, the galaxy sample was split into 
low- and high-redshift sources, and the cosmic shear signal 
was measured on the two subsamples separately. Both sub- 
samples show zero B-modes and the shear amplitude of the 
high-z sample is clearly higher than the low-z one, with 
a ratio in agreement with the cosmic shear predictions. 
The amplitude of the signals from the two subsamples are 
different from each other at all scales with a significance 
level higher than 5-sigma and their shapes follow theo- 
retical expectations of A-CDM dominated universe. This 
strong evidence for the cosmological nature of the signal 
shows that the CFHTLS Deep data will allow us to explore 
the growth rate of cosmic shear signal with redshift, and 



hence the evolution of the dark matter power spectrum as 
function of lookback time. 

Using only Deep data, and marginalizing over h and 
the redshift of sources, we have derived constraints on 
(78 and ilm- We show that the degeneracy between these 
two parameters is partially broken when the analysis is 
combined with data from Wide survey. Assuming flm = 
0.3, we found that as = 0.89 ± 0.06 for P&D and 
CTs = 0.86 ± 0.05 with the halo model, in excellent agree- 
ment with'Van Waerbeke, Mellier fc Ho ekstra 2005| ((78 = 
0.83 ±0.07 ) and Hoekstra et al. 2002 J a>i = 0.86 ±0.05 ). 
Likewise, we derive wq < —0.8 using Deep data alone (see 
IHoekstra et al. 20051 for a deep+wide analysis). 

Our results show that everything is in place to make 
a full scientific use the CFHTLS lensing data, and that 
soon with deeper Deep survey data and wider Wide survey 
data, we will able to provide the best cosmological con- 
straints from weak lensing to date. In particular, we expect 
to explore the growth rate of structure from a tomographic 
cosmic shear measurement, and to better constrain cos- 
mological models from the non-Gaussian features derived 
from a joint analysis of two-point and three-point statis- 
tics. The analysis of three-point statistics in CFHTLS data 
goes beyond the scope of this paper. Indeed, the Deep data 
used for this paper, is not wide enough for such a measure- 
ment. However, three-point statistics will be investigated 
using future samples both for the Deep and Wide survey. 
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